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We study how generic phase transitions taking place in compact stars constructed in the framework 
of the Eddington-inspired Born-lnfeld (EiBf) gravity can lead to anomalous behavior of these stars. 
For the case with first-order phase transitions, compact stars in EiBI gravity with a positive coupling 
parameter k exhibit a finite region with constant pressure, which is absent in general relativity. 
However, for the case with a negative n, an equilibrium stellar configuration cannot be constructed. 
Hence, EiBI gravity seems to impose stricter constraints on the microphysics of stellar matter. 
Besides, in the presence of spatial discontinuities in the sound speed c s due to phase transitions, 

r-f\ | the Ricci scalar is spatially discontinuous and contains 5-function singularities proportional to the 

_ , jump in c\ acquired in the associated phase transition. 
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I. INTRODUCTION 



Based on the gravitational action and nonlinear electrodynamics originally proposed by Eddington [1(, and Born and 

—i ■ Infold [2| , respectively, Bahados and Ferreira recently established a new Eddington-inspired Born-Infeld (EiBI) theory 

£^1 of gravity Q, which reduces to general relativity (GR) in vacuum (see, e.g., 0,[5[ for relevant studies). However, in the 

presence of coupling between matter, EiBI theory demonstrates several distinctive features, including the avoidance 

£~ LV of singularities in the early cosmology and in the Newtonian collapse of pressureless particles, the formation of stable 

i_j. pressureless stars, and the existence of pressureless cold dark matter with a nonzero Jeans length [3l I6l-fl2j . All such 

features stem from the effective repulsive gravitational effect intrinsic to EiBI theory with a positive coupling constant 

OQ , k (see the following discussion for the definition of k) . Hence, in comparison with GR, EiBI theory also allows compact 

£*" ■ stars to attain higher mass limits [g, [ll| . 

Unlike GR, in EiBI theory gravity and matter are nonlinearly coupled together. Therefore, the distinctions between 
GR and EiBI theory arc expected to emerge in the high-density regime. To this end, the static structure and the 
stability of compact stars in EiBI theory have become the focus of several recent papers [g, 0, LL1( ■ As noted in these 
papers, the presence of terms proportional to the spatial derivatives of the density in the equations governing the 
hydrostatic equilibrium and radial oscillation of compact stars in EiBI theory could potentially hamper the related 
studies if nonanalytic (i.e., discontinuous or nondifferentiable) equations of state (EOSs) are considered. To bypass the 
problems arising from such derivatives, in the above-mentioned studies the EOSs employed in constructing com pac t 
stars in EiBI theory are either polytropic or analytic fits mimicking the behavior of realistic nuclear matter [13l |. 
However, even for polytropic EOSs of the form P ex p r , with P, p, and T being the pressure, rest- mass density, and 
a fixed constant, respectively, Pani and Sotiriou have recently discovered that the Ricci scalar diverges at the stellar 
surface if the constant T is greater than 3/2 [14J . Such singular behavior is attributable to the auxiliary field q^ 
proposed in EiBI gravity. In the process of eliminating the auxiliary field in the equation of motion, higher derivatives 
of the matter field emerge and hence lead to a sensitive dependence on the matter distribution. 

In the present paper we investigate how the static structure of compact stars in EiBI theory could be affected by 
drastic changes in the EOS. Specifically, we consider nuclear EOSs with first-order (second-order) phase transitions, 
where the energy density (the sound speed) is a discontinuous function of pressure. Even in GR, the effects of such 
discontinuities are known to be nontrivial [lg, [l|| , and could lead to a third family of compact stars [13, [HI ■ For 
compact stars in EiBI gravity, we find that in the presence of first-order phase transitions, no equilibrium configuration 
exists if the parameter k < 0. Furthermore, the Ricci scalar develops a (^-function singularity and discontinuity at 
the radius where the speed of sound is discontinuous due to phase transitions. We note that such behavior is in stark 
contrast to the case in GR. 

The plan of the paper is as follows. We briefly summarize EiBI gravity and the equilibrium configuration of compact 
stars in Sees. HT1 and IIIII respectively. In Sec. II VI we consider equilibrium stars in EiBI gravity constructed with EOSs 
with first-order phase transitions and investigate the effects and physical implications of such transitions. The singular 
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behavior revealed in the Ricci scalar as a consequence of the discontinuity in the sound speed is studied in Sec. [V] 
Finally, our conclusions are summarized in Sec. I VII We use units where G = c = 1 unless otherwise noted. 

II. EDDINGTON-INSPIRED BORN-INFELD GRAVITY 

Based on the Eddington gravity action and the Born-Infeld nonlinear electrodynamics [l|, [2| , Bahados and Ferreira 
established the EiBI gravity Q, which is summarized here as follows. The starting point of the EiBI theory is the 
action S given by 



r ? 7 = o« QCT ( 9 7<M + d 0°i ~ d ^M) • ( 2 - 4 ) 



S = 16^,/ ^ x [yJ\9^ + KR ^\ 

+Sm[9,^m\. (2.1) 

Here R^ is the symmetric part of the Ricci tensor and is constructed solely from the connection T% . However, the 
matter action Sm depends only on the metric g^ v and the matter field 'I'm- Besides, the determinant of a tensor 
jyn, is denoted by |/ M ^|, and the convention / = |_/" M ^ | is used hereafter. The parameters k and A are related to the 
cosmological constant by A = (A — 1)/k such that, in the limit k — > 0, the action (|2.ip reduces to the Einstein-Hilbert 
action. In the present paper we aim at asymptotic flat solutions and hence consider the case with A = 1. Several 
constraints on the value and the sign of the remaining parameter k have been obtained from solar observations, 
big bang nucleosynthesis, and the existence of neutron stars [3j, |6|, Il2i . Il9| . In particular, for cases with positive re, 
effective gravitational repulsion prevails, leading to the existence of pressureless stars and increases in the mass limits 
of compact stars [g, [ll| ■ 

It is important to note that the spacctime metric g^ v and the connection VZ are treated as independent fields in 
EiBI theory. The field equations are derived by minimizing the action (|2.1j) with respect to g^ and VZ separately, 
and read as Q 

%w = 9V + kR^, (2.2) 

sFk<f v = y/=g<r-Z™y/=g[r' u ' i (2.3) 

with q^y being an auxiliary metric related to the connection, 

1 

2' 
On the other hand, the stress-energy tensor T^ u still satisfies the standard conservation equations, 

V^T^ = 0, (2.5) 

where, as in GR, the covariant derivative refers to the metric g M „. 

It is worthwhile to note that both R^ v and T^ v can lead to a difference in the two metrics g^ v and q^ v . If somehow 
there are some discontinuities in T^ (or its derivatives), then the field equation (|2.3|) implies that g^ v and/or q^ 
will accordingly acquire the corresponding property. However, as R^ w is constructed from the connection Va given 
by (|2.4[) . it contains the second-order derivatives of q^ v . So it is difficult, if not impossible, to maintain the balance 
of the field equation (|2.2[) at these discontinuities if the second-order derivatives of g M „ do not exist. Thus, we expect 
that the discontinuities in T^ v (or its derivatives) will give rise to similar discontinuities in g M „, but not in q^ v . On 
the other hand, in GR the density profile of compact stars (or its spatial derivatives) is in general discontinuous due 
to the presence of phase transitions. Therefore, in the following discussion we consider compact stars in EiBI gravity 
and study specifically how phase transitions could affect the equilibrium configuration and the Ricci scalars derived 
respectively from g^ v and q^. 

III. STATIC EQUILIBRIUM OF COMPACT STARS 

The structure of compact stars in EiBI theory has been studied by Pani et al. [a, 0] and Sham et al. [Il[ . Here we 
briefly review and follow the approach developed in Ref. [Ill ] . For a static and spherically symmetric spacetime, the 
spacetime metric g^ and the auxiliary metric q^ v are taken to be 

g^dx^dx" = -e^dt 2 + e x ^dr 2 + f(r)dn 2 , (3.1) 

q^dx^dx" = -e (3(r) dt 2 + e Q(r W 2 +rW. (3.2) 



The compact star is made of a perfect fluid described by the standard stress-energy tensor 

T» v = (e + PJu'V + Pg ftu , 



(3-3) 



where e and u^ arc the energy density and four-velocity of the fluid, respectively. 

The field equations (|2.2p and (|2.3p lead to a set of relations for the functions </>, A, /, a, and /3 introduced in Eqs. 
(ICT) and d33), 
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(3.5) 

(3.6) 

(3.7) 

(3.8) 



where a = \/l + 8irne and 6 = \/\ — 8ttkP, and hereafter primed quantities denote partial derivatives with respect to 
r. Besides, the conservation of the stress-energy tensor, Eq. (|2.5j) . gives another relation, 



2P' 

P + e' 



(3-9) 



Combining Eqs. (|3.4j) - (|3.9j) . one can obtain two first-order differential equations governing the structure of a compact 
star, 
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the speed of sound c s is calculated from the EOS by c 2 s = dP/de, and the function m(r) is defined by 
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(3.10) 
(3.11) 

(3.12) 
(3.13) 



Equations p.lOj) and (|3.11j> are analogous and in the limit k — > reducible to the well-known Tolman-Oppenheimer- 
Volkov (TOV) equations in GR [2fJ, [21j. With a given EOS P = P(e) and suitable boundary conditions (see below), 
these two equations completely determine the hydrostatic equilibrium configuration of a compact star in EiBI gravity. 
It follows directly from Eqs. (|3.11[) and (|3.13[) that m(r = 0) = 0, and m(r) > increases monotonically with r. 
Besides, it can be shown that 6 > for positive « and 9 ra 8nP{l + 7r«[(e- 3P) 2 /P + 8e] + C(k 2 )} + 2m/r 2 . As the 
typical value of 87r|ft|e at the center of a compact star is less than 0.4 [ll|, for specificity we assume in the present 
paper O > 0. In fact, we have numerically verified that O > except for P -C e. 

The boundary conditions supplementing Eqs. (|3.10[) and p. lip are as follows. First of all, the radius of the star R 
is as usual defined by the condition P(R) = 0. Outside the star where r > R, EiBI gravity is equivalent to GR, and 
g^ v is identical to the Schwarzschild metric (see, e.g., Refs. [22, [23[)- As a consequence, the appropriate boundary 
conditions at the stellar surface r — R arc e = P = 0, a = b = 1, e~ a = e~ A = e 13 — e^ = 1 — 2M/R, where M = m(R) 
is the mass of the star. 

The equilibrium configuration of compact stars has been studied in Refs. [y, l7|, [llj. However, the compact stars 
considered in these references are all characterized by the smooth EOSs P = P(e) where c 2 . = dP/de can be unam- 
biguously defined, and is continuous and nonvanishing except at the stellar surface. Unlike the TOV equations in 
GR, the hydrostatic equilibrium equations (|3.10[) and (|3.1ip demonstrate an explicit dependence on c 2 . Therefore, it 
is reasonable to expect that the behavior of c 2 could strongly affect the equilibrium configuration of a compact star. 
In the following we consider two different cases, namely, (i) c 2 — (as in first-order phase transitions) and (ii) c 2 is 
discontinuous (as in both first- and second-order phase transitions). 



IV. STARS WITH FIRST-ORDER PHASE TRANSITIONS 



For an EOS with a first-order phase transition (see, e.g., Refs. [lg, [24J, |25( ) , there is an interval in which the energy 
density increases while the pressure remains the same, and consequently c 2 = 0. To examine the consequence and 
other related problems of such first-order phase transitions for compact stars in EiBI gravity, we first consider cases 
with c 2 « and expand Eq. (|3.f0[) up to order c 2 , 



dP_ 
dr 



2cVe 



i 



2m 
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(4.1) 



It should be noted that de/dr = cJ 2 dP/dr is still nonvanishing as c 2 — > 0. We study the implications for the cases 
K > and k < separately. 

First of all, for n > 0, we can still get a normal energy density profile in the region with c 2 sa 0, where de/dr is 
approximately equal to a negative constant and dP/dr «0". As shown in Fig.[TJ where a first-order phase transition 
really takes place at a certain pressure (the EOS is adapted from Ref. [2g]), there is a finite region with constant de/dr. 
Therefore the energy density is now a continuous function. Physically, this corresponds to a "mixed phase" region 
(e.g., quark- nucleon mixed phase [lg,[2J]) w ith constant pressure. The thickness of this region is proportional to k, 
and hence vanishes in GR, resulting in a discontinuity in the density profile (see Fig. Q]). The existence of a constant 
pressure shell in compact stars built upon EiBI gravity is attributable to the effective repulsive gravity inherent in 
EiBI gravity with positive k. Actually, the effective repulsion also leads to the formation of pressureless stars, as 
suggested in Refs. [g, 0|- 
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FIG. 1: The density profile of a compact star with an EOS with a first-order phase transition as shown in the inset (P and e 



are in cgs units) is studied in both GR (i.e. k — 0) and EiBI gravity 
density is also fixed at eo 



SvTKeo = 0.1,0.05, with e = 1 X 10 & gcm _d ). The central 



On the other hand, if n < while c 2 m 0, both dP/dr and de/dr will become positive, and in the limit of c 2 — > 0, 
dP/dr — > + and de/dr tends to a nonvanishing positive number. Such behavior is in stark contrast to the case of 
GR where both dP/dr and de/dr are negative definite. If one starts with a certain energy density at the center of a 
star and integrates Eqs. (|3.10[) and (|3.11[) outward, one will find that the energy density drops initially and then rises 
again in the density region where a first-order phase transition occurs and c 2 = 0. As a result, the energy density 
cannot drop to zero and the surface of the star is undefined, implying that equilibrium compact stars characterized 
by EOSs with first-order phase transitions do not exist if k < 0. 

As a matter of fact, for stars with a sufficiently soft phase where c 2 is small, we note a similar problem for EiBI 
gravity with k < 0. We focus on the last term of the modified TOV equation (|3.10p . 



dP 

dr 



P 



(4.2) 



Given an equation of state, this term depends solely on the energy density. When k = 0, i.e., in the case of GR, it 
is always negative. However, when n < and c 2 is small, this term may blow up at certain values of energy density 



regions, and consequently the modified TOV equation cannot be solved. This again implies that for k < 0, compact 
stars with soft EOSs may not be able to support an equilibrium state. 

The predicament encountered in solving Eq. (|3.10j) in the presence of first-order phase transitions can also be 
understood in terms of the concept of the apparent EOS proposed by Delsate and Steinhoff [26[ , who have shown that 
EiBI gravity can be recast as ordinary GR provided that the physical EOS is replaced with an apparent EOS P q (e q ), 
where for perfect fluids the apparent pressure P q and density e q are given respectively by P q = tP + V, e q = re — V, 
with r = [(1 + 87TKe)(l — 87TKP) 3 ] -1 / 2 and V = [r — 1 — kt(3P — c)/2]/k. In first-order phase transitions, where 
dP/de = 0, 

dP^ = STTKJP + e) 

de q a 2 + 3b 2 ' l ' ' 

Therefore, for negative k, dP q /de q < 0. In this case, despite the fact that the TOV equations in GR guarantee that 
P q decreases monotonically towards the stellar surface, the apparent density e q increases instead in the region where 
a first-order phase transition takes place. Hence, by the same token, it is also impossible to construct the equilibrium 
configuration of a compact star governed by the TOV equations and the apparent EOS. 

V. SINGULAR BEHAVIOR OF RICCI SCALAR 

In the above discussion, we see that the equilibrium configuration of compact stars can generally be constructed in 
EiBI gravity with positive k even in the presence of phase transitions. In particular, the energy density e is always a 
continuous function of the radial coordinate r. However, as noted in Ref. jl4| and will be shown in the present paper, 
the associated Ricci scalar still acquires certain singular behaviors (e.g., discontinuities and divergence) due to higher 
derivatives of e. 

In EiBI gravity, the Ricci scalar can be constructed using the metrics q^ and g^, which are respectively given by 

R q = —^TZ [-P' a ' r2 +2/3V + /3'V + 4r/3' - 4ra' - 4e Q + 4] , (5.1) 
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'A' + 20" + 0' 2 + 2(0' " A')y - (^ j J + 4 J- 4 J 



(5.2) 



From the expressions of R q and R g , Eqs. (|5.ip and (|5.2[) . one can show that R q depends on c 2 while R g depends on 
both (? s and (c 2 )' . In the presence of phase transitions of the first or second order, (? s is in general discontinuous at the 
transition point (see, e.g., the inset in Fig. 1). Hence, it is likely that these Ricci scalars could become discontinuous 
or even blow up at radius r t where the phase transition occurs. However, it is readily shown that a, a', /3, j3' and /3" 
are indeed continuous functions of r even in the presence of phase transitions. Hence, R q is still well-behaved even 
in the vicinity of r = ft, which agrees with our intuition as mentioned in Sec. |H] Besides, as shown in Ref. [26| and 
discussed in the last section, R q is indeed the standard Ricci scalar in GR for a star constructed with the apparent 
EOS. Therefore, its regularity is also expected. 

The situation is completely different for R g . It contains a 5- function singularity given by 

R g = ^[Ai(A 2 + 2Aa)]-6(r - n) + ... , (5.3) 

where the notation [F]~ = F(r — rf ) — F(r = rj) for a physical quantity F is introduced, 

A 2 = § + 1, (5.5) 

A 3 = § - ^ (5.6) 

and other finite terms have been suppressed. In addition to the ^-function singularity, R g is also discontinuous there, 
with a jump given by 



[R g ]- = -A ((" V - 2 ^" + oc'P')[A 1 {A 2 + 2A 3 )}_ - P ,2 [A 1 (2A 2 + A 3 )]_ + f3' 2 [A 2 (A 2 + A 2 A 3 + A 2 )]. 
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For reference purpose, the detailed forms of A\ (i = 1, 2, 3) are given below: 

3a 2 + b 2 
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.^ f 8(l + c 2 )c 2 c^ 2\1 
^ kc 2 s y-2A lP '( Z 4-^X (5-9) 






(5.10) 



In general, whereas the 6- function singularity in R g is directly proportional to the discontinuity in c 2 , the jump [R g ]~ 
depends on both [c 2 ]_ and [(c 2 )']_. Thus, for compact stars where first- or second-order phase transitions take place, 
the Ricci scalar R g is singular in the sense that it contains a <5-function-type divergence and spatial discontinuity. As 
is well known, the Ricci scalar in GR is proportional — e + 2>P and is always bounded. The <5-function-type divergence 
in R g discovered here is a distinctive feature of EiBI gravity. 

We have numerically verified the said behavior of R q and R g . For EOSs with phase transitions, either of first- or 
second-order (the EOSs are adapted from Refs. [lg,|25J|, respectively), R q is found to be a continuous function, while 
R g (see Fig. [2]) is discontinuous and indeed blows up at the place where there is a jump in c 2 . The discontinuity in 
R g obtained numerically in fact agrees with the expression given by Eq. (|5.7|) . 

VI. DISCUSSION 

In this paper, extending our previous work [11], we have studied compact stars in EiBI gravity using realistic EOS 
models with phase transitions and discovered several anomalies in the behavior of such stars. For EOSs with first- 
order phase transitions, compact stars in EiBI gravity with k > exhibit an anomalous "mixed phase" region where 
c 2 = 0, the pressure is a constant and the energy density e is still a continuous function of the radius r. This is in 
contrast to the situation in GR where e is discontinuous at the transition point. It is the effective pressure inherent in 
EiBI gravity with positive k that leads to the existence of the constant pressure region. On the other hand, for k < 0, 
the equilibrium configuration for compact stars with first-order phase transitions (or soft enough EOSs) cannot be 
constructed because e never vanishes and hence it is not possible to fulfill the boundary conditions (e = P = 0) at 
the stellar surface in order to match the interior solution smoothly to the Schwarzschild spacetime. 

As a side remark, in general for k^O, we find that one cannot construct the simplest quark star model, described 
by the MIT bag model EOS, since e is finite when P = for this EOS. Note that this does not pose a problem in GR 
because one only requires P — at the surface in GR. 

In GR, one can always construct a compact star model by solving the TOV equation with a prescribed EOS. The 
EOS model is required to satisfy only a few general conditions which are already imposed by the microscopic theory 
(e.g., the stability condition dP/de > and the causality limit c s < c). On the contrary, EiBI gravity appears to 
put a more severe restriction on what kind of EOS one can use to build a stellar model. It is true that whether 
compact stars in nature exhibit a phase transition in their interiors is still a matter of debate. Nevertheless, it seems 
unreasonable that the underlying gravitational theory would put a constraint on which microscopic EOS model one 
can use to construct a theoretical stellar model in the first place. Even in the situation where a compact stars with a 
first-order phase transition can be constructed in EiBI theory (i.e., when k > 0), we still need to face the S- function 
singularity of the Ricci scalar R g as we have shown in Sec. [V] 

Our discovery reported here complements the recent work of Pani and Sotiriou [lj] in which they demonstrated 
the singularity of R g at the surface of a polytropic sphere for any T > 3/2. In fact, the singular behaviors of R g 
discussed in the present paper and Ref. [lj] stem from a common physical origin, namely the emergence of the first- 
and second-order derivatives of the pressure in the Ricci scalar R g . These derivatives in the matter field arise as a 
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FIG. 2: The Ricci scalar R g (solid line, left-scale) and energy density e (dashed line, right-scale) of a compact star where a 



first- (upper panel, central density = f .2eo ) or second-order (lower panel, central density 
plotted against the radius r. 87TKeo is equal to 0.1 for both cases 
lines in the figures) emerge at the places where c a is discontinuous. 



3eo) phase transition occurs are 



For R g , S- function singularities (as indicated by the vertical 



consequence of the elimination of the auxiliary metric q^ defined in Eq. (|2.4p . (Sec Ref. [TJ for a detailed discussion. 
) They are absent in the GR case and lead to enigmatic singularities in the Ricci scalar, as mentioned above. 

In conclusion, EiBI gravity is appealing because it can avoid some of the singularities that plague GR by introducing 
nonlinear coupling between matter and gravity. However, it is also the same nonlinear coupling that leads to the 
anomalies of compact stars with phase transitions, as we have found in this paper. This renders the viability of EiBI 
gravity questionable. 
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